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ABSTRACT 



Context. Determining ages of the outer and the inner lobes of so-called double-double radio galaxies (DDRGs) is crucial 
for understanding the active cycles of galactic nuclei, the phases of interruption of the jet flow, and physical conditions 
in the surrounding galactic and intergalactic medium governing the jets' propagation. A recognition and understanding 
of these conditions during the restarted jet activity is of special interest. 

Aims. We determine the ages and other physical characteristics of the outer and the inner lobes of the DDRG 
J1548— 3216, as well as the properties of the surrounding environment during the original and the restarted phase 
of the jets' activity. 

Methods. Using the new low-frequency and high-frequency radio images of this galaxy, we determined the shape of the 
spectrum along its lobes and performed the classical spectral-ageing analysis. On the other hand, we applied the ana- 
lytical model of the jet's dynamics, which allowed us to derive the physical conditions for the source's evolution during 
the original jet propagation through the unperturbed IGM, as well as those when the restarted new jet propagates 
inside the outer cocoon formed by the old jet material that passed through the jet terminal shock. 
Results. The dynamical age estimate of the outer and the inner lobes is 132±28 Myr and ~9±4 Myr, respectively. The 
synchrotron age in the outer lobes systematically rises from ~25 Myr in the vicinity of the lobes' edges to about 65-75 
Myr in the centre of the old cocoon. These ages imply an average expansion speed along the jets' axis: (0.012±0.003)c 
in the outer lobes and (0.058±0.025)c in the inner lobes, but the latter speed would be ~0.25c when they were of age 
less than 1 Myr. We find that the jet power during the restarted activity is about ten-fold fainter than that of the 
original jet. Similar disproportion is found for the internal pressures and the magnetic field strengths in the old cocoon 
and those in the inner lobes. This disproportion can be effectively reduced by assuming the same equations of state for 
the emitting particles and the magnetic fields within the old and the new lobes. However, we think that our assumption 
of the non-relativistic equation of state for the old cocoon and the relativistic one for the new lobes is more justified. 

Key words, galaxies: active - galaxies: evolution - radio continuum: galaxies 



^ ■ 1. Introduction 



Although the intrinsic time evolution of powerful radio 
sources of Fanaroff-Rilcy type II (FRII; Fanaroff & Riley 
1974) is largely understood and described with a num- 
ber of analytical models (e.g. Kaiser, Dennett-Thorpe 
& Alexander 1997; Blundell, Rawlings & Willott 1999; 
Manolakou & Kirk 2002; Kino & Kawakatu 2005), there are 
still several unanswered questions about the duty-period of 
the active galactic nucleus (AGN), the jet production pro- 
cesses, its interaction with the external gaseous environ- 
ment including the intergalactic medium (IGM), and the 
contents of the radio lobes as a part of the low-density 'co- 
coon'. 

The double-double radio galaxies (DDRGs) are char- 
acterized by two pairs of unequally-sized edge-brightened 
(FRII-type) lobes sharing the same radio core. In most 
of them the outer and inner double structures are aligned 
well. The existence of such radio sources is the evidence 
that the jet activity in AGN may be not continuous during 
the lifetime of a source. In fact, an intermittent produc- 
tion of jets can be connected with stochastic transitions 



between two accretion modes: the standard one - with an- 
gular momentum transmitted outwards by viscous torques 
within the accretion disk - and the 'magnetic' one, with 
the developed large-scale magnetic fields and related MHD 
winds (Nipoti, Blundell & Binney 2005; Kording, Jester 
& Fender 2006). Sikora, Stawarz & Lasota (2007) incor- 
porated the above idea into the spin paradigm scenario. 
Postulating that the efficient production of relativistic jet 
requires both a large black hole (BH) spin (as in the model 
of Blandford & Znajek 1977) and an efficient collimation 
mechanism (cf. Begelman & Li 1994), they noted that the 
intermittent jet activity observed in active galaxies accret- 
ing at high rates may be due to intermittent collimation of 
the central Poynting flux-dominated (so called 'Blandford 
& Znajek') outflow by heavier and slower MHD wind gen- 
erated in the inner parts of the accretion disk undergoing 
state transitions. In the framework of this interpretation, 
the jet axis in the subsequent jet activity epochs is expected 
to be the same, since this axis is determined by the spin of 
the central BH, which should not change substantially on 
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short (<C 100 Myr) timescales (see the discussion in Sikora 
et al. 2007). 

Such an interrupted production of jets is evidently im- 
printed in the radio morphology of DDRGs. We are inter- 
ested in certain aspects of these sources. Are the ages and 
internal densities of the inner lobes much lower than these 
values for most 'normal' radio sources of similar physical 
size as the inner lobes? Is the density of the pre-existing co- 
coon much lower than in the unperturbed galactic and IGM 
environment, or is this density higher? The first case would 
strongly suggest that the new inner structure is formed in 
a channel drilled through the old cocoon by the former jet 
activity cycle, which has been modelled by the numerical 
MHD simulations of Clarke & Burns (1991). They predict 
that 'the restarted jet will always be overdense (denser than 
its immediate surroundings) if the original jet is under- 
dense relative to the quiescent IGM'. While the restart- 
ing jet model accounts for many of the observations, there 
remain some profound discrepancies difficult to be recon- 
ciled (cf. Clarke 1997). The second case would imply an 
efficient replacement of the inner lobes by the heavier ex- 
ternal medium (e.g. Kaiser, Schocnmakers & Rottgcring 
2000; Brocksopp et al. 2007). 

The radio galaxy J1548-3216 (PKS B1545-321) is a re- 
markable example of DDRG in which the newly restarted 
jets propagate through the remnant cocoon of a previ- 
ous active phase (Subrahmanyan, Saripalli & Hunstead 
1996; Saripalli, Subrahmanyan & Shankar 2003, hereafter 
referred to as SSS2003). This galaxy has recently been ex- 
tensively studied by Safouris et al. (2008, hereafter referred 
to as S2008), especially under the aspect of an observa- 
tional constraint for the 2D and 3D numerical simulations 
of the restarted jet provided by Clarke & Burns (1991) and 
Clarke (1997), respectively. In S2008, the authors suggest 
that observational data are consistent with a picture that 
the restarted jets generate narrow-bow shocks, and the in- 
ner lobes in this galaxy are a mixture of cocoon plasma 
re-accelerated at the bow shock and new jet material re- 
accelerated at the termination shock. They propose that 
the evolution of the restarted jets and the inner lobes is 
strongly influenced by an entrainment of the external IGM 
into the pre-existing cocoon. 

In this paper the spectral-ageing and dynamical analy- 
ses of J1548— 3216 are performed with the aim of (i) deter- 
mining the synchrotron age distribution in the outer lobes 
of this galaxy (in the old cocoon), (ii) estimating of this 
age in the inner lobes, (iii) comparing these ages with the 
dynamical ages estimated with the DYNAGE algorithm 
of Machalski et al. (2007), (iv) determining the jet pow- 
ers during the first and the second phase of activity, as 
well as other physical parameters characterizing the lobes 
and their environments, such as the particle density, en- 
ergy density, internal pressure, magnetic field strength and 
its density, propagation speeds of the lobes along the jets' 
axis, etc. Most working approaches in this paper are sim- 
ilar to those applied in our previous publications on other 
giant-sized radio galaxies: DDRG J1453+3308 (Konar et 
al. 2006) and a further ten selected galaxies (Jamrozy et 
al. 2008; Machalski, Jamrozy & Saikia 2009, hereafter re- 
ferred to as MJS2009). The analyses presented in this paper 
are based on the new radio observations recently conducted 
with the Giant Metrewave Radio Telescope (GMRT) and 
Very Large Array (VLA), and on the Australia Telescope 
Compact Array (ATCA) and VLA archival data kindly pro- 



vided to us by Vicky Safouris and Ravi Subrahmanyan. The 
new observations and the data reduction are presented in 
Sect. 2. The resulting total-intensity 334, 619, 1384, 2495, 
and 4860 MHz total-intensity images are used in Sect. 3 
to derive radio maps of the outer double structure of the 
investigated galaxy, as well as to extract the inner dou- 
ble structure from a background of the underlying cocoon 
formed during the earlier phase of the nuclear activity. The 
spectral-ageing analysis of the outer and the inner struc- 
tures is described in Sect. 4. while the dynamical analysis 
is presented in Sect. 5. Results of these analyses, as well 
as our contribution to the aspects of the restarted nuclear 
activity, the environmental conditions ruling the new jets 
propagation within the relict cocoon, and their energetics 
- derived with another approach than applied in the pre- 
vious studies of this radio galaxy by SSS2003 and S2008 - 
are discussed in Sect. 6. 

For the purpose of calculating the linear size, volume, 
and luminosity of the lobes, we use cosmological parameters 
ft m =0.27, fi vac =0.73, and H Q =71 kms^Mpc- 1 . 

2. Observational data and their reduction 

The observing log for all the observations is listed in 
Table 1, which is arranged as follows. Columns 1 and 2 
show the name of the telescope and the array configuration 
for the former and the recent VLA observations; columns 3 
and 4 give the frequency of observations and the primary 
beamwidth; columns 5 and 6 show a typical angular resolu- 
tion and an rms noise level achieved in the resulting images 
of the radio galaxy investigated. The last column 7 gives 
the dates of the observations. More details of these obser- 
vations are given below. 

The low-frequency GMRT observations at 334 and 619 
MHz were made in the standard manner, with each obser- 
vation of the target source interspersed with observations 
of calibrator sources. The phase calibrators B1714— 252 (at 
334 MHz) and B1626-298 (at 619 MHz) were observed 
after each of several 20 min-lasting exposures of the tar- 
get centred on the core position. 3C286 was used as the 
flux density and bandpass calibrator based on the scale of 
Baars et al. (1977). At each of the two frequencies the to- 
tal observing time on the target source was only about 150 
min because of very limited observing time scheduled for 
the project. Unfortunately, a large part of 334-MHz data 
were strongly affected by radio frequency interference, and 
these data had to be flagged in the reduction process, which 
further reduced the quality of the data. Acceptable data 
were edited and reduced with the NRAO AIPS package. 
All these data were self-calibrated to produce the best pos- 
sible images. 

At the frequencies of 1384 and 2495 MHz, the archival 
data taken with the ATCA and VLA arrays are used. 
In particular, the ATCA 2495 MHz map of the total 
structure published by SSS2003, as well as the combined 
ATCA+VLA 1384 MHz of the total structure and the high- 
resolution VLA 1384 and 4910 MHz map of the inner dou- 
ble published by S2008. For the purpose of specifying a 
high-frequency spectrum of diffuse lobes of the outer double 
structure, we made other 4860 MHz observations of the tar- 
get source with the VLA in its DnC configuration. Again, 
3C286 and B1522— 275 were used for the amplitude and the 
phase calibrations, respectively. Two 20 min exposures of 
the fields centred on each of the two outer lobes were re- 
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Table 1. Observing log. 
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Teles- 


Array 


Obs. 


Primary 


Ang. 


rms 


Observing 


cope 




freq. 


beam 


resol. 


noise 


date 






[MHz] 


[arc min] 


[arc sec] 


[mjy beam ] 




GMRT 




334 


80 


15 


0.24 


2008, May 24 


GMRT 




619 


43 


7.5 


0.12 


2008, Mar 8 


ATCA+VLA* 




1384 


40 


3.5 


0.04 


2001, Dec 2 


VLA* 


BnA 


1384 


30 


3.5 


0.04 


2002, May 31 


ATCA* 




2495 


24 


5 


0.05 


1999, Mar 7 


VLA 


DnC 


4860 


9 


15 


0.03 


2008, Jun 5 


VLA* 


BnA 


4910 


9 


1.3 


0.02 


2002, Jun 2 



*: Archival data; courtesy of Vicky Safouris and Ravi Subrahmanyan. 



duccd, self-calibrated, and combined into one image of the 
entire source. 

3. Observational results 

3.1. New radio images 

A full-resolution GMRT 334 and 619 MHz images, as well 
as the VLA/DnC 4860 MHz image, are presented in Fig. 1. 
Our new images, especially those at low frequencies, con- 
firm the overall morphology of J1548— 3216 already pre- 
sented and discussed by SSS2003 and S2008, i.e. that both 
the outer lobes are edge-brightened and rather sharply 
bounded. Likewise in those papers, our images also do not 
show any evident hot spots or very compact structures at 
the ends of the lobes, and both low-frequency images con- 
firm a distinct pair of emission peaks along a bright rim 
at the western end of the NW lobe. As in the archival 
ATCA and ATCA+VLA data, the inner double structure 
(a pair of relatively bright, narrow lobes) is strongly im- 
mersed into the diffuse bridge of emission extended from 
the bright edges of the outer lobes towards the radio core. 
A flare of the bridge transverse to the source's axis in the 
vicinity of the core, very well shown at the ATCA 1384 
MHz image in SSS2003, is also pronounced at both low 
frequencies. This flare is missing at the 4860 MHz image, 
which suggests a very steep radio spectrum in that part of 
the structure (cf. Sect. 3.4). However, a missing flux den- 
sity at this frequency is negligible. Indeed, the area marked 
with the dashed line in the right panel of Fig. 1 is about 60 
restoring beams, and a missing flux is likely between the 
rms noise level and the first contour CI in this image, both 
multiplied by ^60 beams, i.e. between ~1.8 mjy and ^6 
mjy. Even the missing flux of 6 mjy will be about 1.3% 
of the total flux density of 449 mjy given in Table 2. Such 
a loss docs not affect the spectral analysis performed in 
Sect. 3.4. 

However, the new radio images, showing nothing espe- 
cially new in respect to the archival ones, extend the ob- 
servational data from two to the five different frequencies 
ranging from 334 MHz to 4860 MHz. This is the neces- 
sary and sufficient condition for performing the undertaken 
analyses. 

3.2. Extraction of the inner double structure 

To perform the spectral-ageing and dynamical analyses sep- 
arately for the outer and the inner structures, we have to 
extract emission of the inner lobes from the underlying 
'background' radiation of the outer lobes. At the observing 



frequencies of 334 and 619 MHz this is made by exclud- 
ing of the visibility data taken with baselines shorter than 
2 kA and 3 kA, respectively, while at 4860 MHz visibili- 
ties with spatial frequencies less than 2.5 kA are excluded. 
This effectively resolved out a large part of the underly- 
ing bridge's emission. Somewhat different approach was ap- 
plied at 2495 MHz. Having the final ATCA image at this 
frequency only but not the original UV data, we used the 
AIPS task IM2UV which allows a Fourier transformation 
of the image reconverting the data back to a UV data file. 
Then a similar procedure, as described above, was applied 
to the reconverted 2495-MHz UV data excluding the visi- 
bilities at baselines shorter than 2.5 kA. 

The resulting 334, 619, 2495, and 4860 MHz images 
of the inner double structure are shown in Fig. 2. The 
corresponding archival ATCA 1384 MHz and VLA/BnA 
4910 MHz images are included for comparison. Besides the 
brightest parts with the leading heads of the inner lobes, 
our new images confirm the presence of another weak emis- 
sion region in the inner NW lobe detected by SSS2003 in 
their 2495-MHz total intensity images of the inner struc- 
ture. Unfortunately, the dynamic range of our images is too 
low to detect more of the connecting emission seen in their 
Fig. 6. We estimate that such a missing flux is from about 
10% at 334 MHz to about 2% at 4860 MHz. All these im- 
ages, except VLA/BnA, brought to a common scale using 
the AIPS task HGEOM and convolved to the angular res- 
olution of 16" x 16", are used to make a longitudinal section 
along the inner structure. Such 'slices' at the five observing 
frequencies are shown in Fig. 3. To avoid problems with any 
missing flux, we restrict our spectral and dynamical analy- 
ses of the inner lobes to their brightest regions indicated in 
Fig. 3. 

3.3. Outer lobes cleaned from the inner structure 

To analyse physical properties of pure outer lobes of 
J1548— 3216, the inner double structure was subtracted 
from the images of the entire radio source partly shown 
in Fig. 1. For this purpose, all of those images were also 
brought to a common scale (a map size, cell size, coor- 
dinates of the map centre) and convolved to the angular 
resolution of 16" x 16". The images of the inner structure 
were blanked over regions outside the extracted inner lobes 
and then subtracted from the convolved maps of the en- 
tire source using the AIPS task COMB. The net images of 
the outer lobes (rotated by 35°) are shown in Fig. 4. On 
the first of these images, a division of the radio structure 
into 18 strips, each of them 28" wide, is shown. The first 
plotted contour on the ATCA+VLA image is exceptionally 
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Fig. 1. Full-resolution GMRT 334 and 619 MHz, as well as VLA DnC-array 4860 MHz images of the entire structure of 
the radio galaxy J1548— 3216 (PKS 1545-321). The first contour level, CI, is given in each image. The contour levels are 
(1, 2, 4, 8,...)x CI mJy/beam. In all the images, the restoring beam is indicated by an ellipse. The cross indicates the 
position of the optical parent galaxy. The dashed outermost contour on the VLA image encloses an area of the flare of 
the emission bridge within which a missing 4.86 GHz flux density is estimated in the text. 
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J1548-3216: inner double 
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Fig. 3. Longitudinal section along the inner double struc- 
ture. The horizontal brackets indicate regions of the struc- 
ture subject to the spectral and dynamical analyses. 



high (about 8x rms noise level) to clear it from spurious 
jagged contours that appeared after the convolution of the 
original map with the beam of 16" x 16". The integrated 
flux densities measured in the consecutive strips and plot- 
ted vs distance of the strip's centre from the core position 
(the strips' centres are separated by the angular distance of 
28"/ cos 35°) form a longitudinal section along the cleaned 
outer structure shown in Fig. 5a. The flares or spurs in the 



central region of the outer structure increase the total flux 
density in the strips S7, S8 and Nil, N12 causing its peaks 
marked Fl, F2, and F3 in Fig. 5a. The brightness peaks of 
the leading heads of the new inner lobes lying at the ra- 
dio axis (indicated with the dotted vertical lines), almost 
coincide with the positions of these strips. A spectral steep- 
ening and the spectral age within the strips are analysed in 
Sect. 4.1. 



3.4. Radio spectra 

The integrated flux densities of the total source, as well as of 
its outer and inner lobes, are given in Table 2. All columns 
are self-explanatory; outNW and outSE indicate the NW 
and SE lobes of the outer double structure, while innNW 
and innSE - the NW and SE inner lobes, respectively. 
Because a spectral fit, especially with the SYNAGE, is very 
sensitive to a lack of low-frequency data (cf. MJS2009), the 
160-MHz flux densities of the outer lobes are estimated by 
subtracting 300 mJy (assumed flux density of the inner dou- 
ble at this frequency based at a spectral index of about 0.6, 
cf. Sect. 4.2) from the total flux density measured with the 
Culgoora array, and dividing the net flux between the two 
lobes in a proportion similar to those observed at the higher 
frequencies. 

Distributions of the low-frequency af|| 4 and high- 
frequency «4860 spectral index vs distance from the core 
measured along the axis of the outer structure cleaned from 
the inner lobes are shown in Fig. 5b. The wavy ridge and its 
side flares do not show any peculiarity in the spectral index 
distribution shown in Fig. 5b. Both the low-frequency and 
the high-frequency indices exhibit a systematic steepening 
from the heads of the outer lobes towards the centre. The 
low-frequency spectral index rises from ~0.5 to ~0.9, while 
the high-frequency one steepens from ~0.8 to ~2.0 at the 
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Fig. 2. GMRT, ATCA and VLA images of the inner double structure of J1548— 3216. The first contour level is given in 
each image. The contour levels are (1, 1.41, 2, 2.83, 4, 5.66,...)xCl mJy/beam. The restoring beam is indicated by an 
ellipse. The cross marks the position of the optical galaxy. 



evident depression of emission at the centre of the bridge. 
Such a large continuous steepening of the spectra suggested 
a systematic increase in the synchrotron age of relativistic 
particles enclosed in the old cocoon, i.e. an increase from 
the lobes' head towards their flaring ends. 

We do not attempt to analyse a distribution of the syn- 
chrotron age in directions transverse to the main axis of 
the source, hence spectral index distributions over the en- 
tire area of the outer structure are beyond the scope of this 
paper. Nevertheless our data show similar spectral features 
as those seen in the map of the spectral index aflgl in 
SSS2003 (their Fig. 3), i.e. the steepest spectra appear at 
eastern ends of the strips S7 and S8, and at western ends 
of the strips Nil and N12. We cannot confirm a distinctly 
steeper spectrum along the southwestern edge of the outSE 
lobe appearing in their map, but at least something similar 
is not pronounced in the spectral index cv^f 4 or agfg 4 . 



3.5. The radio core 

The J2000.0 position of the radio core determined from the 
high-resolution images is R.A.: 15 h 48 m 58 s . 05 and Dec: 
— 32°16'56".9, which is less than 1 arcsec away from cen- 
tre of the parent galaxy imaged with the Anglo- Australian 
Telescope (AAT) by S2008. The flux densities of the core, 
measured on the images presented in this paper, are col- 
lected in Table 3. These flux densities suggest a mildly in- 
verted spectrum without a sign of time variability. 



Table 3. Flux densities of the radio core 



Freq. 


Score 


Observing 


[MHz] 


[mJy] 


date 


619 


1.52 


2008, Mar 8 


1384 


2.17 


2002, May 31 


2495 


2.50 


1999, Mar 7 


4860 


2.53 


2008, Jun 5 


4910 


2.57 


2002, Jul 30 



4. Spectral ageing analysis 

Remembering all the serious problems with both the prin- 
ciples and the practical application of spectral-ageing cal- 
culations to physical conditions in radio sources described 
in detail in MJS2009, the spectral age in different parts of 
the lobes, i.e. the time elapsed since the radiating particles 
were last accelerated, is determined using the classical the- 
ory that describes the time evolution of emission spectrum 
of a single population of particles with an initial power-law 
energy distribution (e.g. Myers & Spanglcr, 1985; Carilli et 
al., 1991). The initial energy distribution of the relativistic 
particles is a power-law function, n^^dji = n '-f~ p dji, of 
their initial Lorentz factor, 7$. The power p corresponds to 
the initial (injection) spectral index a; n j, which can be, in 
principle, estimated from the observational data until the 
synchrotron frequency of the minimum electron Lorentz 
factor lies far outside the observable low-frequency spec- 
trum. Fortunately, a spectral turnover at low frequencies 
is not observed in the radio spectra of the extended FRII- 
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Fig. 4. Radio images of the outer lobes cleaned from the inner double structure and convolved to the angular resolution 
of 16" x 16". The first contour level is given in each image. The contour levels are as in Fig. 2. The division into 18 
28"-wide strips, used for the spectral steepening and spectral age analysis, is shown in the first of these images. 

type radio sources. On the other hand, the spectral break chrotron) age, T syn , and the magnetic field strength, B, 
frequency above which the radio spectrum steepens from through 
the injected power law, ^br, is related to the spectral (syn- 

r syn [Myr] - 50.3-^^-K(l + *)]" 1/2 , (1) 
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Table 2. Flux densities of the total structure (Total) and the outer (outNW and outSE) and inner (innNW and innSE) 
lobes of J1548-3216. 



Freq. 


Total 


outNW 


outSE 


innNW 


innSE 


Ref. 


[MHz] 


[mjy] 


[mjy] 


[mjy] 


[mjy] 


[mJy] 




(1) 


(2) 


(3) 

/onnnxicn\* 


(4) 


(5) 


(6) 


(7) 



334 
619 
843 
1384 
1400 
2495 
4860 
4910 



4926±740 
3274±265 
2519±200 
1815 ±90 
1842±764 
1019 ±31 
449 ±25 



2095±314 
1383±111 

820±41 

479±15 
198±20 



2642±396 
1758±141 

913±46 

484±15 
217±22 



58.0±9.2 
43.4±4.6 



130.0±19.7 
87.6±7.6 



26.5±3.3 53.3±4.0 



17.5±3.0 
11.8±1.6 
11.9±0.9 



35.5±3.2 
19.7±2.2 
23.1±1.5 



(5) 
(5) 
(2) 
(4) 
(3) 
(4) 
(5) 
(4) 



References: (1) Slee 1995; (2) Jones & McAdam 1992; (3) NVSS (Condon et al. 1998); (4) Flux densities measured on the 
archival ATCA and VLA images (Saripalli et al. 2003; Safouris et al. (2008); (5) this paper; (*) estimated division of the total 

160 MHz flux density between the outer lobes, cf. the text. 
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J1548-3216: outer structure 
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Fig. 5. a) Integrated flux densities in the consecutive strips 
along the outer double structure vs distance from the core. 
The vertical lines indicate positions of brightness peaks in 
the inner lobes, b) Low-frequency and high-frequency spec- 
tral indices in these strips. 



the observed radio spectra using the SYNAGE algorithm 
of Murgia (1996). 

4.1. The outer structure 

4.1.1. Determination of a; n j and i^r values 

To determine the value of ai n j, we fit the CI, CIE, and 
JP models to the flux densities of the entire outer lobes 
(given in columns 3 and 4 of Table 2) treating a; n j as a free 
parameter, and realizing that decidedly the best fit to the 
data is achieved with the JP model. Fits of the JP model of 
radiative losses to the flux densities of the outer SE and NW 
lobes are shown in Fig. 6. The values of ai n j=0.583lQ;Q7Q 



and ainj =0.540 



af|g| indices of 



09 ^ found by the fit correspond to the 



0.05 

^0.7 and ^0.6 previously determined by 
SSS2003 for the brightest regions at the SE and NW heads 
of the outer structure, respectively. These fitted «i n j indices 
are used to determine values of Vhr in the 18 parallel strips 
covering the entire outer structure of the radio source. The 
JP models of the spectra within these 18 strips are collected 
in Fig. 7. A distance of the strip's centre from the core, the 
resulting value of ^br, and the relevant value of Xred giving a 
goodness of the fit in each of 18 strips, are given in columns 
2, 3, and 4 of Table 4, respectively. 

4.1.2. Determination of magnetic field strength values and 
the spectral ages 

In consistency with the approach applied in our previous 
spectral-ageing analyses of giant radio galaxies (Jamrozy 
et al. 2005, 2008; Konar et al. 2006), the magnetic field in 
Eq. (1) is identified with an 'equipartition field', B cqv , which 
provides equipartition between the total energy densities 
of the relativistic particles and the magnetic field (u c « 
ub). The required values of B cqv are computed with Miley's 
(1980) prescription for the general formula 



Bcqv oc (1 + k) 



V 



2/7 



(2) 



where i?ic=0.318(l±z) 2 is the magnetic field strength 
equivalent to the inverse-Compton microwave background 
radiation. B and Bic are expressed in units of nT, while 
i^br is in GHz. Values of a- m j and ^ r are found by the fit to 



where k is the ratio of the energy content of relativistic pro- 
tons to that of electrons (adopted as k—1), L is the luminos- 
ity of a given strip calculated by integration of its spectrum 
from a frequency equivalent to a minimum Lorentz factor, 
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Fig. 6. Spectra of the outer lobes fitted with the JP model, as described in the text. 



7min ~ 1 for the relativistic electrons to the upper limit of 
100 GHz, and V is the volume corresponding to that slice. 
The derived values of B cqv and the resulting spectral ages, 
T syn , are given in columns 5 and 6 of Table 4, respectively. 

The distribution of this spectral age vs distance from 
the core measured along the axis of the outer structure is 
shown in Fig. 8. 



4.2. The inner structure 

The spectrum of each of the two inner lobes, i.e. the flux 
densities given in columns 5 and 6 of Table 2, is fitted with 
the CI model. The fits (shown in Fig. 9) suggest a similar 
initial slope of the spectrum of both the lobes of ~ 0.6±0.1 
and the spectral break of 161^^ 1 ' ) GHz and 44^2 n ^ GHz 
for the NW and SE lobes, respectively. (The formal ±lcr 
errors are enormous due to the practically straight spec- 
tra.) The volume of the lobes is calculated assuming their 
cylindrical geometry with a minimum angular size 24" x 7" 
(height x base diameter) for the NW lobe and 42" x 8" 
for the SE lobe where these dimensions are measured in 
the VLA image of S2008 (their Fig. 5). In this case, the 
equipartition magnetic field strength, calculated with the 
prescription of Milcy (1980), is 0.63±0.16 nT and 0.65±14 
nT, respectively. Using these values, a 'mean' spectral age of 
the radiating particles in the lobes is 5.4+4 g Myr for innNW 
lobe and lO.l+g^ Myr for innSE lobe. However, adopting 
the full length of the lobes as the cylinder's height, i.e. 
74. "6 and 89. "5 for innNW and innSE lobes, the magnetic 
field strengths reduce to 0.46 ± 0.10 nT and 0.43 ± 0.10 nT, 
while the ages increase to 7.0±l° 5 Myr and 13.9±ll 3 M y r > 
respectively. The resulting spectral ages of the inner lobes 
are discussed in Sect. 6.1. 



5. The dynamical age analysis 

This analysis is performed using the DYNAGE algorithm 
of Machalski et al. (2007) which is based on the analytical 
model for the evolution of FR II type radio sources, combin- 
ing the pure dynamical model of Kaiser & Alexander (1997) 
with the model for expected radio emission from a source 
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Fig. 8. Spectral age of relativistic particles in the outer 
structure cleaned from the inner lobes plotted vs distance 
from the radio core. 



under the influence of the energy loss processes published 
by Kaiser et al. (1997, known as the KDA model). With 
this algorithm we derive the dynamical age of the lobes t, 
both the outer and the inner ones, the effective injection 
spectral index a- m - } , which approximates the initial electron 
continuum averaged over a very broad energy range and 
over the present age of the source, the jet power Qj et , and 
the central density near the radio core po, which determines 
the ambient density in which the jet propagates. 

A detailed description of how to apply the above algo- 
rithm is published in MJS2009. It is worth explaining here 
that determining of values of these four free parameters of 
the model is possible by a fit to the observational parame- 
ters of a source (or its lobes): its projected linear size D, the 
volume V, the radio power P v and the radio spectrum a u , 
which provides (P v )i a t a number of observing frequencies 
z=l, 2, 3,... As in the KDA model, we assume a cylindri- 
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are given in Table 4. 



cal geometry of the lobes (cocoon), thus V = 7rZ? 3 /4i?^ 
where Rt is their axial ratio. The values of the few other 
free parameters of the model have to be assumed. These 
are the central core radius cto, the exponent (3 describing 
the ambient density profile in the simplified King's (1972) 
model p(d) = po(d/ao)~" , the Lorentz factors determining 
the energy range of the relativistic particles used in integra- 
tion of their initial power-law distribution 7i, m i n and 7i, m ax, 
the adiabatic indices of the three 'fluids' with individual en- 
ergy densities: the jet material, rj et , the magnetic field, Tb, 
and the ambient medium, T x (cf. Kaiser et al. 1997). Since 
the emitting region consists of these three fluids, the model 
also takes the adiabatic index of the lobe (cocoon) into ac- 
count as a whole, T c . The two other free parameters we 
have to assume are k' - the ratio of the energy density of 
thermal particles to that of the relativistic electrons - and 
8, the orientation of the jet axis to the observer's line of 
sight. Following KDA, in the DYNAGE algorithm the as- 
sumed energy equipartition is expressed by the ratio of the 
energy densities of the magnetic field and of the particles, 



C = (l+p)/4 = (1- 
whole source are ao= 



ai n j)/2. The values adopted for the 



10 kpc, f3=1.5, 7;, 



=1, 7i, 



=10' 



and 6>=90°. A decrease in 6 to -70° (cf. S2008) would result 
in ~ 6% increase of D and ~ 7% increase of t (cf. Eq. (6) in 
Section 5.3). As we are interested in an age difference be- 
tween the lobes rather than in their absolute age value, the 
latter one is less important. The values of T c , Tb, T x , and 
k' assumed for the outer and the inner lobes are given in 
next sections. The observational data of these lobes, used in 
the DYNAGE fitting procedure, are given in Table 5. Most 
columns are selfexplanatory, the entries in columns (6)-(9) 
give the ratios of the size and luminosity of the given lobes. 



Given the values of a- m j, Q 



inj, Vjeti P0 



and t, several other 



physical parameters of the source can be specified, e.g. the 
internal pressure in the lobes p c (t), their energy density 
u c (t), a ratio of the kinetic energy delivered by the jet to 
the energy radiated out, (Qj c t x t)/(u c x V), and an average 
expansion speed of the lobes, D/(cxt). The assumption of 
the energy equipartition condition allows the magnetic field 
density Ub(£) and the field's strength B(t) to be estimated. 
The detailed expressions were given in our previous paper 
(cf. MJS2009). The table with all notations for physical pa- 
rameters used through the paper is given in the Appendix. 
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Table 4. Break frequency, equipartition magnetic field strength, and spectral age of emitting particles in consecutive 
strips through the outer lobes of J1548— 3216 (cf. Fig. 4). 
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However, the age and other physical parameters, fitted 
independently for either lobe of a given double source, may 
be significantly different; i.e., any difference between the 
fitted values of a parameter is found to be greater than 
the uncertainty of the fits. This is a consequence of the 
usual asymmetries between the lobes in their length and 
luminosity. The difference arises if the same density pro- 
file of the ambient medium along the opposite lobes is as- 
sumed. The ratios between these parameters of the lobes of 
J1548— 3216, both the outer and the inner ones, are shown 
in columns (6)-(9) in Table 5. On the other hand, we can 
expect that Qj ct and p have the same values in the so- 
lutions for the opposite lobes, since they characterize an 
energy-emitting process in the central AGN. Also a large 
difference in age is rather unlikely. Therefore, following a 
similar ageing analysis in MJS2009, we consider the inde- 
pendent solutions, as well as the two self-consistent solu- 
tions for the opposite lobes, hereafter denoted as solutions 
A and B. 



5.1. Independent solutions for the individual outer lobes 

The images in Fig. 4 suggest that the diffused outer lobes 
may comprise a fraction of thermal particles, thus the lobes, 
as a whole, may have a non-relativistic equation of state. 
Therefore, we assume r c =rB=r x =5/3 and A;'=10. The lat- 
ter value is less than ~ (25-140) calculated by S2008 with 
the assumption that the hot spot pressure, ph s , equals a 
minimum pressure p m - m , though they expected that a true 
hot spot pressure should be higher, i.e. phs > Pmin- Indeed, 
when studying the large-scale X-ray environment of se- 
lected FRII radio sources, Belsole et al. (2007) find that the 
internal pressure in their lobes, p c , is ~(l-5) times higher 
that the minimum (equipartition) pressure. Considering 
that the ratio Phs/Pc in the DYNAGE algorithm varies from 



^4 to ^20 (depending on the value of Rt), the assumed 
value k' is justified. 

The model solutions, i.e. the parameter values resulting 
from the independent fits, are listed in columns 2 and 3 of 
Table 6. 



5.2. Self-consistent solutions for the outer lobes 

The differences between the values of the model parameters 
for the opposite lobes found in the independent solutions 
come from different environmental conditions (and/or dif- 
ferent magnetic fields) on both sides of the core region. For 
this reason, in the first kind self-consistent solution (so- 
lution A), we averaged the values of Qj e t,out and po.out 
found for the opposite outer lobes (given in columns 2 
and 3 of Table 6), and now treat them as the fixed free 
parameters of the model, (<3j t,out}=T-115 x 10 38 W and 
(poout)=5.445 x 10~ 22 kgm -3 , respectively. Given these val- 
ues, we can determine another value of (5 for each of the 
two lobes, hereafter denoted as /3 SC .A- To do that (follow- 
ing MJS2009), we equalize values of the ambient density at 
the head of the outer lobes resulting from the independent 
solution and the self-consistent solution A: 

Pa,out = Po,out(D/a o )~ « (po,out)(D/a o )~ 0scA , hence 

log((po,out)/Pa,out) 



P. 



Sc.A 



Iog(D/oo) 



(3) 



Given a value of (3 S c.a, the expected age of a lobe in the 
frame of the self-consistent solution is 



^out ~ — 



(Qjct.out) 



(4) 



J. Machalski et al.: Spectral and dynamical ages of GRG J1548— 3216 



11 



T 1 1 I I I I I 




Frequency [MHz] 



Fig. 9. Spectra of the inner lobes fitted with the CI model, as described in the text. 

Table 5. Observational parameters of the outer and the inner lobes used to fit the dynamical model. The monochromatic 
powers are expressed in WHz _1 sr~ 1 ). 
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2.96 


3.37 
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1.20 


Rt 


2.7 


3.1 


11.8 


10.6 










log P334 


24.756 


24.654 


23.470 


23.115 


19.32 


34.59 


1.26 


2.26 


log P 6 l 9 


24.585 


24.481 


23.300 


22.991 


19.27 


30.90 


1.27 


2.04 


log Pi384 


24.310 


24.263 


23.085 


22.780 


16.79 


30.41 


1.11 


2.02 


log P 48 60 


23.716 


23.668 


22.705 


22.432 


11.46 


17.22 


1.12 


1.68 



where c\ is a dimensionless constant dependent on the val- 
ues of (3, Tj, T x , and T c , given by Eq. (25) in Kaiser & 
Alexander (1997). The model parameter values resulting 
from the self-consistent solution A are listed in columns 4 
and 5 of Table 6. 

The data in Table 6 show that the difference between 
the lobes' ages inferred from the solution A is greater than 
the found in the independent solution. This is not what 
we would expect for the actual ages of the opposite lobefl 
Another alternative, self-consistent solution is plausible in 
which these ages are very similar (especially if we suspect 
that an orientation of the jets' axis in giant radio galax- 
ies is close to 9 ~ 90°), and any differences between the 
linear extent and luminosity of the lobes come from an in- 
homogeneity either in density distribution of the ambient 
gaseous environment or in magnetic field. Since significant 
differences between the jet power and the radio core param- 
eters in the opposite directions along the jets' axis are not 
plausible, in the self-consistent solution B we assume the 
same values of (Qjot.out) and (po.out) for both outer lobes 
(as in solution A) and the same age (i out } = 132 Myr, i.e. a 
mean of the ages determined in the independent solution, 
and ao=10 kpc. In such a scenario, a value of (3 sc .b can be 



1 In MJS2009 we showed that for some sources the solution 
A can diminish the age difference. This usually happens if the 
shorter lobe is brighter than the larger one 



calculated from Eq. (4) substituting (i ut) for t out and (3 sc .b 
for Psc.a- As a result, 

Pscb(D) = (5) 

This solution does not give an unequivocal result for the fit, 
because ci is a rising function of (3. As shown and discussed 
in MJS2009, the parameter space within which opposite 
lobes of a significantly asymmetric source would have the 
same age is usually large. Nevertheless, an example of a 
reliable solution B for (5 sc .b (calculated with Ci(f3—1.5)), 
the corresponding value of a m j , and other model parameters 
of the outer SE and NW lobes are shown in columns 6 and 7 
of Table 6, respectively, for a comparison. It is worth noting 
that in this solution the product (Qjct.out) x (^out) does not 
provide a minimum of the jet kinetic energy that is assured 
in the independent solutions. However, such a minimum is 
always shallow (cf. Machalski et al. 2007) and the above 
objection is not very important. 

The data in Table 6 show that the values of the model 
free parameters do not differ much in the three solutions 
considered. However, the identity of Qj C t,out, A), out, and t ou t 
in the opposite outer lobes postulated in the solution B 
result in a greater difference between the values of p a .out 
and j? c ,out than in the remaining solutions. In particular, 
it suggests more than twice denser ambient environment 
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around the head of the outer NW lobe than around the 
head of the opposite SE lobe. These physical conditions 
seem to be supported by the presence of a distinct pair of 
emission peaks along the bright rim at the end of the NW 
lobe, while a similar emission is absent in the SE lobe. Also 
the mean pressure, p c ,out in the NW lobe is about 30% 
higher than in the SE lobe, and their ratio found in the 
solution B is the highest. We discuss this point again in 
Sect. 6.2. 

5.3. Independent solution for the inner lobes 

In the case of the inner structure, we assume that (i) the ob- 
served emission arises from the narrow lobes (cocoon) , not 
the restarted jets; (ii) the jets' and lobes' material has a rel- 
ativistic equation of state with rj et = Tb = T c = 4/3 with 
no thermal particles, thus fc'=0; and (hi) the restarted jets 
propagate within rarefied and uniform (with /3=0) medium 
of the relict outer cocoon formed by the old jets' mate- 
rial that passed through the jet terminal shock. Since the 
observed spectra of the inner lobes show no curvature be- 
low the frequency of 4.9 GHz, especially for the innNW 
lobe where the SYNAGE fit suggests above 20 GHz (cf. 
Fig. 9), the DYNAGE algorithm will not be able to find 
a unique solution for the dynamical age, i.e. to determine 
values of Qj C t,inn, /Oo,inn, and t lnn , even if a value of a- m j is 
known. Its formal fit with the SYNAGE is 0.606^;^ and 

0.611±°;°^ for the innSE and innNW lobes, respectively. 
Therefore for the purpose of DYNAGE calculations, we as- 
sume here that a maximum value of a- ln j cannot exceed the 
values of 0.606 and 0.611, but can be as low as 0.510 and 
0.525, respectively. Moreover, the width of the inner cocoon 
can be larger than the lobes' widths determined from the 
images in Fig. 2, therefore we admit a twice larger width for 
these lobes corresponding to i?x=6.5 (instead of 11.8 and 
10.6) supposing that the best age solution for the innSE 
lobe lies within the model space parameters limited from 
one side by the values a; n j =0.606 and i?T=H-8, and from 
the other side by ai n j =0.510 and Rt—Q-5. For the opposite 
innNW lobe, the limiting pairs of the model parameters are 
cv inj =0.611, i? T =10.6 and cv inj =0.525, i? T =6.5. 

The sets of solutions resulting from the fit of the model's 
free parameters to the linear size and the radio powers of 
the inner lobes (given in columns (4) and (5) of Table 5), 
are presented in Fig. 10. This diagram clearly shows that 
the spaces of model parameters for the opposite inner lobes 
do not overlap. Obviously the lobes' asymmetries in the 
luminosity and size are too large to allow a comparable age 
and jet power solution in the model. Moreover, a selection 
of adequate pair of Qj e t,inn and t lnn values is not possible 
until a value of po,inn is fixed by means of some additional 
constraint. 

Let us therefore consider the limiting values for the 
core density /9o,inn within the old outer lobes. On the one 
hand, the upper limit for the cold gas density may there- 
fore be provided by studies of the internal depolarization 
of radio emission produced by the extended lobes of FRII- 
type radio galaxies. For example, Garrington & Conway 
(1991) found that the product of the cold gas number den- 
sity and the lobes' magnetic field strength is on average 
n g x B < 0.5 x 10 3 m" 3 nT. This, with the B oqv « 0.64 nT 
determined for the inner lobes in Sect. 4.2, gives roughly 
Po.inn ~ WpTig < 1.3 x 10~ 24 kgm~ 3 . We note in this con- 



text that the above equipartition magnetic field strength is 
compatible with the typical values 0.3 nT< B < 3 nT found 
by means of multiwavelength analysis of the non-thermal 
lobes' emission (e.g. Kataoka & Stawarz 2005; Croston et al. 
2005). However, the ambient gas density within the old co- 
coon of DDRGs is likely lower than that of the typical FRII- 
type sources with linear sizes comparable to those charac- 
teristic for the inner double structures. It can be supposed 
that p ,inn < Pa,out, where p a , out = (po,out)(A>ut/ao) _:L5 - 
With (p ,out) ~ 5.4 x 10~ 23 kgm~ 3 and D out w 500 kpc, we 
have po,mn < 1.6 x 10~ 25 kgm~ 3 . This value, corresponding 
exactly to the baryon density of 10 -4 cm~ 3 (10 2 m~ 3 ) typ- 
ical for the intracluster medium (ICM) of nearby (z < 0.2) 
clusters of galaxies (Croston et al. 2008), can be considered 
as a very upper limit for po.inn- 

On the other hand, the lower limit is provided by an 
amount of matter injected into the outer lobes by the old 
jet, which is 

^ Qjct,out^out 



where m p , n p , and n c are proton mass, proton, and electron 
number densities, respectively, and k = n c /n p . Allowing 
for mildly relativistic value of Vj ct on large scales and only 
slight electron dominance in terms of number density within 
the outflow (Sikora & Madejski, 2000; Celotti & Ghisellini, 
2008), one can expect roughly rj et (l + k) < 10, and thus 
the lower limit for the gas density within the old cocoon 
Po,mn > 10~ 30 kgm~ 3 . However, such a density is extremely 
low and insufficient to explain the observed characteristics 
of the inner structure. Indeed, power of the jet propagating 
through the medium with po.inn i<5 10 -30 kgm~ 3 and f3 = 
(calculated from suitably transformed Eq. (4)) sufficient to 
expand a lobe with Rt « 11 {c\ « 4) to the length D « 150 
kpc with the speed Vh c, i.e. during t <; 0.5 Myr, is 
~ 5 x 10 35 W. Therefore, in the framework of the KDA 
model, the power radiated out at the frequency of 1.4 GHz 
will only be about ~ 1.4 x 10 18 WHz _1 sr _1 , i.e. almost five 
orders of magnitude lower than the observed one! Kaiser has 
(2000) already argued that some 'additional material must 
have passed from the environment of the source through 
the bow shock surrounding the outer source structure into 
the cocoon'. 

Perhaps a better lower limit for the gas density re- 
sults from the numerical simulation of an intermittent jet 
by Clarke & Burns (1991). They show that the restarted 
jet must propagate through a thcrmalized material left by 
the old jet, and this medium is about 40 times less dense 
than the initial ambient medium. In our case po,inn > 
Po,out(-Dinn/ao) _1 ' 5 /40. Therefore we estimate that 2.5 x 
10~ 26 < po.inn < 1.6xl0~ 25 kgm~ 3 . Taking a mean of these 
extreme values, (po,mn) = 6.3 x 10~ 26 kgm~ 3 , we achieve 
further limits for the space of parameters in Fig. 10. The 
vertices of the solid-line- marked diamonds in Fig. 10 indi- 
cate Solutions Of iinn and Qjet.inn for (po,inn) = 6.3 x 10~ 26 

kgm~ 3 and the four combinations between the values of 
a; n j and Rt. (i) 0.611 and 10.6; (ii) 0.611 and 6.5; (iii) 
0.525 and 10.6, and (iv) 0.525 and 6.5 for the innNW lobe 
for instance. In particular, the marked areas indicate that 
the dynamical age of the inner lobes is within a range of 
6-13 Myr, which is in accordance with the synchrotron age 
derived in Sect. 4.2. 
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Table 6. Fitted physical parameters of the outer lobes, with brackets ( ) showing the values assumed within the given 
solution, (cf. the text). 



Parameter Indepen. solution Self-consist. solut.A Self-consist, solut. B 

outSE outNW outSE outNW outSE outNW 



(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


P, {P SC .A)(l3sc. B ) 




1.50 


1.50 


1.494 


1.507 


1.586 


1.402 






0.552 


0.531 


0.536 


0.545 


0.527 


0.554 


Qjct,o Ut (xl0 38 W) 




1.211 


1.019 


(1.115) 


(1.115) 


(1.115) 


(1.115) 


po,out(xl(T 23 kgHr 


■ 3 ) 


5.585 


5.297 


(5.445) 


(5.445) 


(5.445) 


(5.445) 


Pa,out(xl0 -25 kgm" 


- 3 ) 


1.521 


1.570 


1.535 


1.574 


1.040 


2.382 


tout (Myr) 


142 


122 


146 


119 


(132) 


(132) 


Pc,out(xl(T 13 Pa) 




1.472 


1.459 


1.393 


1.552 


1.291 


1.687 


Bout(nT) 




0.322 


0.320 


0.313 


0.330 


0.300 


0.345 


Uout =u c x V(xl0 53 J) 


1.884 


1.200 


1.784 


1.280 


1.624 


1.399 


D/(t out x c) 




0.0118 


0.0130 


0.0115 


0.0133 


0.0128 


0.0120 



The value of (/Oo,inn) adopted in our calculations is 
more than two orders higher than its estimate in S2008 
(< 5 x 10 -28 kgm -3 ) and close to the ionized gas density 
in clusters of galaxies (cf. Croston ct al. 2008). This seems 
to contradict the finding of SSS2003 that the investigated 
radio galaxy is not in a cluster, and one could expect that 
the gas density inside the old cocoon is much lower than 
the value (/Oo.inn) adopted in the above calculations unless 
an efficient entrainment process occurs. However, an order 
of magnitude decrease in this value, i.e. to about 6 x 10 -27 
kgm -3 , will result in a two-fold decrease in the estimated 
age of the inner structure only, i.e. to about ^3 Myr. 

Nevertheless, there is a way to avoid any assumption 
about the environment's density. It is not an easy task to 
detect and measure the high-frequency radio spectrum of 
the inner lobes. The calculation of the flux density expected 



from the model with 



'inj 



0.56, 



t jet, inn 



= 1.26 x 10 37 W, 



A),hm = 6.8 x 10 -2b kgm -3 , R T = 11, and t = 9 Myr are 
20.6 mJy, 11.0 mJy, 4.6 mJy, and 1.9 mJy at v = 10.6, 30, 
90, and 230 GHz, respectively. Supplementing the observed 
spectrum with the above flux densities and fitting the CI 
model of energy losses, we find i^r = 110.4 GHz, which 



with B cqv — 0.45 nT gives r f 



syn 



8.6 Myr. 



Similar to the outer lobes, the independent solutions for 
the inner lobes cannot provide comparable jet powers for 
the opposite lobes. The diamond-shaped areas in Fig. 10 
are significantly separated along the abscissa. We therefore 
consider an alternative solution of their age. 



5.4. Alternative solution 

An explanation of the difference in Qj C t,inn values by differ- 
ent external density profiles along the opposite new jets is 
problematic if we assume a constant density with (3=0. The 
nearer opposite sides of the 'diamonds' correspond to the 
largest acceptable differences between the parameters Rt 
and cvi n j . Even completely different values of these parame- 
ters for the opposite inner lobes do not allow even a partial 
superposition of these diamonds in the (log) plane ti nn - 
Qjct,inn- However, as pointed out and discussed in MJS2009, 
a plausible equalization of Qj ot values can be achieved by 
changing the ratio between the energy densities of magnetic 
field and of particles, For example, the six- fold decrease 
of ( in the innNW lobe results in desired increase in Qj t,inn 
in this lobe. This alternative solution is marked in Fig. 10. 
We discuss this solution in Sect. 6.2. 




36.0 36.5 37.0 37.5 

10g(Q Je t.in„[W]) 



38.0 



Fig. 10. Model solution for the inner lobes. The two pairs of 
dotted diagonal lines for each of the lobes indicate the age 
vs jet-power relation corresponding to the limiting values 
of the lobe's axial ratio, Rt, and injection spectral index, 
cVjnj, marked in the diagram. The diamond-shaped solid- 
line areas show credible dynamical age solutions for the 
two lobes. The dashed diamond indicates an alternative 
solution for the innNW lobe providing an equalization of 
the jet powers for the opposite lobes (cf. the text). The 
solid vertical line indicates the value of Q- Jct adopted for 
the alternative age solution. 



Adopting the jet power, (Qj e t,inn) = 1-26 x 10 37 W, 
which intersects the diamond-shaped areas in the middle in- 
dicating the independent solution range for the innSE lobe 
and the alternative solution for the innNW lobe (Fig. 10), 
and considering that (po.inn) = 6.3 x 10 -26 kgm -3 , as well 
as /3=0, a probable age of the inner lobes estimate, ti nn , is 
given by Eq. (4) where the parameters denoted with 'out' 
are replaced by those denoted with 'inn'. Given two differ- 
ent values of cl corresponding to two extreme values of Rt 
adopted here for either lobe, the age estimate is between 6.3 
and 8.8 Myr, and between 7.9 and 12.0 Myr for the innNW 
lobe and innSE lobe, respectively. 

A number of physical parameters fitted with the model 
for the inner lobes of J1548— 3216, with the ages as above, 
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are provided in Table 7. Columns (2) and (3) give their 
values for the innSE lobe with i?T=H-8 and innNW lobe 
with i?T=10.6, respectively, while columns (4) and (5) give 
their values for these lobes but with i?x=6.5. Although the 
fitted values of Qjct.inn and poinn differ a little from the val- 
ues (Qjet.inn) - 1.26 x 10 37 W and (p ,inn) = 6.3 x 1CT 26 
kgm -3 used to estimate the age from Eq. (4), their ratio 
Poinn/Qjct,inn always equals (po,inn}/(Qjet,inn). This occurs 
because of the property of the DYNAGE algorithm allow- 
ing determination of the values of both p and Qj et by the 
fit of the model free parameters to the observed size and 
power of a given source (lobe). The last two lines in Table 7 
give the expansion velocity of the jet's head resulting from 
the model and exhibit its deceleration with time. Since the 
length of the jet is 



it) 



Cl 



Qjet 



1/(5-/3) 



f 



3/(5-/3) 



implying that (6) 



dr 3 ci 
«(*) - 37 = 



dt 5-/3 \p a , 

1/3 / \ (/3-2)/3 



\ V(5-/3) 

if! t (/3-2)/(5-/J) ) or(7) 



v(r) 



3ci Q 



jet 



5-/3 



(8) 



The next to last line in Table 7 gives the expansion velocity 
at the assumed radius of the radio core do, while the last 
one gives this velocity at the actual length of the given inner 

lobe Ann- 



6. Discussion of the results 

In the sections below, we discuss the derived ages and phys- 
ical parameters of the source, both its outer and inner dou- 
ble structures, as well as those associated with the original 
unperturbed ambient environment and the relict cocoon. 



6.1. Spectral age of the outer and the inner structures 

As expected, the synchrotron age of both the outer lobes, 
T syn , increases with distance from their edges. It starts from 
~25 Myr and rises to about 65-75 Myr in the centre of the 
old cocoon. The lowest estimated age, far from the zero 
value, can be related to the time that passed from the 
last acceleration of the emitting particles at the region of 
the old jet interaction with the original unperturbed IGM. 
However, in a number of similar spectral-ageing analyses 
(cf. Liu, Pooley, Riley 1992; Jamrozy et al. 2005, 2008) 
a mean of the ratio between the youngest and the old- 
est emitting particles is always about 0.40 ± 0.15 and is 
weakly dependent on the size and/or the age of sources. 
It means that the youngest synchrotron age is found usu- 
ally between 5 Myr and 20 Myr for FRII-type sources with 
D £ 100 kpc and t £ 20 Myr (cf. Machalski, Chyzy & 
Jamrozy 2004). Similar values (~ 18.7 Myr and ~ 20.3 
Myr) have been found in the detailed ageing analysis of the 
DDRG J1453+3308 (Konar et al. 2006), for which the time 
that passed since the jets stopped actively feeding the hot 
spots was estimated by Kaiser et al. (2000) as about 1.1 
Myr. We return to this point again in Sect. 6.3. 



The evidently curved age-distance tracks in Fig. 8, es- 
pecially for the outNW lobe, support the expected decel- 
eration of the jet head's advance speed and probable back- 
flow present in giant radio galaxies (e.g. Schoenmakers et 
al. 1998; Lara et al. 2000). Indeed, the advance speed of 
the jet's head given by Eq. (8), v(r) = «h s , and a backflow 
velocity, v D f, will form opposite vectors. If a quotient of 
separation between the strips in Table 4 and a difference 
between the derived spectral (synchrotron) ages in these 
strips, Ar/Ar syn = (v S yn), is more than v(r) at a given dis- 
tance from the core, a backflow is present. If (v D f) = e(t>h s ), 
then (v syn ) = (1 + e)(vh s ). For example, the entries in 
Table 4 show that, in outSE lobe at r = 332 kpc, e w 2.2 
and Vbi ~ 0.024 c, while closer to the core at r = 131 kpc: 
e w 1.3 and v hi w 0.017 c. In outNW lobe at r = 300 kpc, 
e w 2.0 and v hf « 0.024 c, while at r = 98 kpc: e « 0.4 and 
Ubf « 0.005 c. 

A spectral age for the new inner lobes is quite prob- 
lematic. Since this structure is relatively smaller and much 
fainter than the outer structure, we were not able to cut it 
into strips and to check whether and how much the spec- 
tral index varies along the lobes. The spectral index map 
derived from the the ATCA+VLA and VLA images at 1384 
and 2495 MHz, published by S2008 (their Fig. 5), does not 
indicate any steepening that would signal a spectral age- 
ing. Though the flux densities measured in the entire inner 
lobes (cf. columns (5) and (6) in Table 2) allowed the CI 
model fit to the data and an estimate of spectral breaks in 
the observed spectra, the fit's uncertainty is so large that 
the resulting 'mean' spectral ages within the lobes are very 
uncertain as being formally in a range from ~0.1 Myr to 
~14 Myr. Taking an uncertainty of the Rt value of the in- 
ner lobes into account (which determines their volume, V , 
used to estimate B eqv ), and that even the axial expansion 
of the lobes with the speed of light would require the time 
of ^0.5 Myr, we may adopt that a spectral age of these 
lobes is about 9 ± 5 Myr. 

6.2. Dynamical age of the outer and the inner structures 

Following similar ageing analysis made by MJS2009 for a 
limited sample of ten giant radio galaxies, we solved the 
dynamical age and other physical parameters of the two 
outer lobes independently and explored the resulting for- 
mal differences between those fitted parameters that should 
actually not differ in the opposite lobes, like the jet power 
and the central core density. In this way we can perform 
a deeper search for properties of the surrounding ambient 
medium and conditions of the jets' propagation during the 
active phases of the jet production in the common central 
engine. In the two kinds of proposed self-consistent solu- 
tions, we require identical values of Qj ct and po for the op- 
posite lobes both of the outer and of the inner structures. 
However, in MJS2009 we showed that the fits with such 
an assumption results in either decreasing or enlarging the 
age difference between the lobes, depending on the actual 
asymmetries in the lobes' length and luminosity. The self- 
consistent solution A for the outer lobes of J1548— 3216 en- 
larges the age difference from ~20 Myr in the independent 
solution to ~27 Myr (cf. the entries in columns (2)-(5) in 
Table 6). As the above differences are too large to be caused 
by a projection of the source on the sky, we consider the 
solution B in which we demand the same age for both lobes. 
Taking a mean of the age values in the independent solu- 
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Table 7. Fitted physical parameters of the inner lobes for different values of the lobe's axial ratio Rt- 
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Parameter 


Indepen. solution 


Alternat. 


solution 




innSE 





innNW 


/ -i \ 

(i) 


(2) 


(3) 


(4) 


(5) 




11.8 


6.5 


10.6 


6.5 


Oinj 


0.510 


0.605 


U.02o 


O.DlO 


Qjet,inn(xlO' :> ' W) 


1.53 


1.21 


1.39 


1.34 


„ / w i n — 26 1 „ . — 3 \ 

Poinn(,xlO kgm ) 


7.70 


6.07 


7.17 


6.13 


tinn(Myr) 


7.9 


12.0 


6.3 


8.8 


Pc, l nn(xlO- 14 Pa) 


2.91 


3.34 


3.45 


4.10 


B inn (nT) 


0.20 


0.22 


0.22 


0.24 


C/inn = U C X \/(xl0 5 ° J) 


1.50 


5.66 


1.27 


3.44 


D/(i x c) 


0.072 


0.047 


0.075 


0.053 


i>(r = ao)/c 


0.286 


0.189 


0.267 


0.189 


V(r = Ann)/C 


0.043 


0.028 


0.045 


0.032 



tion, i.e. (i ut)=T32 Myr, we find it as satisfactory for the 
outer lobes. This is worth emphasizing that this age fully 
agrees with its estimate in S2008, (30-200) Myr, though 
they derive it from another consideration about the energy 
budget. 

As expected, the self-consistent solution B suggests dif- 
ferent exponents (3 sc .b in the power-law density profiles and 
different effective injection spectral index ai n j, where [3 sc .b 
value (according to Eq. (6)) governs the lobe's length, while 
ctinj value influences the energy losses and the resulting 
brightness of a given lobe. Different fitted values of (3 sc .b 
and a; n j for the opposite outer lobes result in about twice 
higher IGM density at the edge of the outer NW lobe as 
compared to that at the SE lobe. As a result, the bright rim 
and the two 'warm' spots in the outNW lobe may indicate a 
higher pressure ratio between the lobe's head and the IGM, 
in other words, this lobe may be more overpressured than 
the opposite outSE lobe. 

The age for the inner structure is estimated by assuming 
a scenario in which the new restarted jets propagate within 
rarefied and uniform (Ann=0) medium of the old outer co- 
coon formed by the old jet material that passed through 
the jet terminal shock and eventually partly mixed with 
the surrounding thermal gas (cf. the discussion in Kaiser 
et al., 2000). The assumption of uniformity is justified by 
the fact that the sound crossing time within the outer lobes 
is much longer than the dynamical timcscales of evolution 
of the outer structure. In the case of the inner structure, 
the age solution is formally undefined due to almost pure 
power-law radio spectrum of the new lobes. On the one 
hand, such a lack of the spectrum steepening can tell us 
that the structure is rather young; on the other hand, it 
prevents determination of its quantitative value (cf. also 
SSS2003). In the framework of the DYNAGE algorithm the 
inability to specify an explicit age solution means that the 
space of the model's free parameters (<5j e t, Po> t, ^inj) is 
very large (the space between two diagonal lines for a given 
lobe in Fig. 10). Within this space, kinetic energy of the jet, 
Qjet x t, is almost constant, so that the criterion of its min- 
imum (a crucial aspect of the algorithm) is not applicable. 
As shown in Sect. 5.3, the above space can be reduced a 
little by fixing values of some of these parameters, e.g. the 
value of po (hence the diamond-shaped areas in Fig. 10). 

However, the straight, unaged spectra of the inner lobes 
are not just a single factor precluding age determination. 
The other one is asymmetry both in their length and bright- 
ness, which causes the spaces of parameters allowed for the 



opposite lobes to not overlap each other. The assumption 
of uniform density of the pre-existing cocoon with /3—0 pre- 
cludes a differentiation of this parameter for the lobes, thus 
we pointed out another parameter that also rules the age 
solution, namely the ratio between energy densities of the 
magnetic field and of the particles, £. The calculations in- 
dicate that a significant departure from the cquipartition 
conditions in one of the inner lobes (but not in the opposite 
one) enables a relative equalization of the jet powers and 
attainment of comparable ages for both lobes. Therefore, 
a combination of the alternative solution for the innNW 
lobe and the independent one for the innSE lobe is required 
to attain comparable values of the three parameters in- 
volved: QjoUnn, Po, inn, and i; nn . Although this alternative 
age solution for the innNW lobe is only a numerical result 
of the fitting algorithm, such a departure seems to be sup- 
ported by X-ray observations, which allow measurements of 
electron energies of radio lobes and magnetic fields. In the 
lobes of the well-studied nearby radio galaxies Centaurus B 
and Fornax A, the inferred electron energy densities exceed 
those of magnetic fields by a factor of 5 - 8 (cf. Tashiro et 
al. 1998, 2009; Isobe et al. 2006). 

Of course, the above combined age solution is odd, in 
which a significant departure from equipartition conditions 
happens in only one lobe of the inner double structure 
evolving in supposedly uniform medium, poann ~ const. 
Perhaps, these different conditions might arise when the 
magnetic field within the old cocoon is strongly inhomoge- 
neous giving rise to enhanced observable radiation in the 
high- field regions (cf. Eilek et al. 1997; Kaiser 2005). An in- 
homogeneity of the field may indicate a faint blob of emis- 
sion detected at all the observed frequencies between the 
radio core and the bright end of the innNW lobe (cf. Fig. 2). 
Understanding an expected criticism of the above specula- 
tions, we argue that the estimated dynamical age of the 
inner lobes of the investigated radio galaxy, ti nn w 9 ± 4 
Myr, is plausible. This age accords with its spectral age 
estimate given in Sect. 5.3. 

6.3. Kinetic energy and internal pressure 

Kinetic energy delivered by the original jets to the old co- 
coon during the time of ~132 Myr is found to be ~10 54 J, 
while the energy radiated out is ~ 1.3 x 10 53 J (cf. entries in 
Table 6). The latter value is about twice higher than esti- 
mated in S2008 using another approach to the energy bud- 
get. One of the important results of our dynamical analysis 
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is the power of the old and the new jets. We find that the jet 
power during the previous phase of the nuclear activity was 
~ 1.2 x 10 38 W, while the power of the new jets (although 
we could only determine a range of its value ~ 5 x 10 36 W 
to ~ 3 x 10 37 W; cf. Fig. 10) seems to be almost one order 
less than the former one, but in a good agreement with its 
range of [1.1, 4]xl0 44 ergs" 1 if 55° < 6 < 78° estimated by 
S2008 (cf. their Fig. 16). Such a high ratio Qj e t,out/Qjet,inn 
impairs results of the ageing analysis of a few other DDRGs 
published by Kaiser et al. (2000) in which they assumed the 
same Qj e t values during two successive phases of activity. 
The difference between the outer and inner jet power found 
in our analysis may indicate that the spin of the central 
BH cannot be the only factor determining the total power 
of the relativistic outflow, since this spin is not expected to 
change significantly in a relatively short timescale, which 
would be involved. However, we note that, in the original 
Blandford & Znajek (1977) model, the power of the outflow 
extracted from the ergosphere of the rotating BH depends 
not only on the BH spin, but also on the magnetic field 
intensity within the innermost parts of the accretion disk, 
which in turn should depend on the accretion rate (or, pos- 
sibly, on the accretion mode). Thus, the difference between 
the power of the original and the restarted jet in the radio 
galaxy J1545— 3216 may still be reconciled with the modi- 
fied spin paradigm (Sikora et al. 2007; Sikora 2009) if the 
accretion disk undergoes significant structural changes be- 
tween different activity epochs in a single source. 

The DYNAGE algorithm clearly shows that a demand 
of Qjot.out ~ Qjct.inn would lead to an unlikely solution of 
the model for the inner structure with the lobes' age U nn « 
0.8 Myr, the old cocoon density po,inn ~ 3.3 x 10~ 28 kgm~ 3 , 
and the average expansion speed Di nn /t lnn « 0.61c! One 
can conclude that until an independent measure of po.inn is 
available (e.g. from internal depolarization of radio emission 
and/or X-ray inverse-Compton emission from the region of 
the inner lobes), a numerical algorithm like the DYNAGE 
will not be able to specify more accurate values of Qj e t,inn, 

^inn, -^inn/^inn, etc. 

The internal (cocoon) pressures, p c , resemble the jets' 
disproportion. The pressure in the old cocoon, derived from 
the fit, is about five times higher than the pressure estimate 
within the new younger lobes (cf. the entries in Tables 5 
and 7). However, this is a result of the different equations 
of state assumed for the jet, particle, and magnetic field 
'fluids': non-relativistic for the old cocoon and the relativis- 
tic one for the inner double structure. If (3=0, the cocoon 
pressure is 

PcK (^) 2/5/M ~ 4/5 - 

For a lobe with given D and J?t we find from Eq. (6) that 
(Qjet/po) xi 3 = const[m 5 ]. Thus 

p c oc p t~ 2 . 

Assuming the non-relativistic equation of state for a lobe 
at age t, the model requires a higher jet power and pre- 
dicts proportionally the higher core density needed to fit 
the lobe's length and power, higher than the values of Qj et 
and p fitted in the relativistic case. As a result, in the 
lobe at age t, p c oc po- The calculations show that even a 
large decrease in the lobe's age gives little increase of p c , 



so that the pressure achieved in the non-relativistic case is 
not attainable. 

6.4. Magnetic fields 

With the standard assumption about the energy equiparti- 
tion between the relativistic particles and magnetic fields, 
the derived strength of uniform magnetic field in the outer 
lobes, -B cqv , slightly varies from ^0.33 nT at the edges of the 
lobes to ~0.24 nT in the centre of the old cocoon. Similar 
values are derived from the magnetic energy density, uq , via 
the dynamical considerations. For the inner lobes, the field 
estimates are incompatible with the values derived with the 
Miley (1980) formula during the spectral-ageing analysis, 
which are almost twice higher than the corresponding val- 
ues implied by the dynamical analysis. Our calculations (in 
the framework of the KDA model) show that, to increase 
strength of magnetic field in the inner lobes to its level de- 
rived from the spectral analysis, it is necessary to assume 
the non-relativistic equations of state for the jet's material 
and the magnetic field 'fluids', as well as to admit some 
admixture of thermal particles. Obviously such an assump- 
tion would level out the propagation conditions during the 
initial and the restarted phases of activity, which is rather 
unacceptable. 

7. Conclusions 

The new low-frequency and high-frequency radio contin- 
uum maps of the double-double radio galaxy J1548— 3216 
(PKS 1545—321) are used to complement its already pub- 
lished maps at 22 and 13 cm (Saripalli, Subrahmanyan & 
Shankar 2003; Safouris et al. 2008) and to perform both 
the spectral-ageing and the dynamical analysis of this re- 
markable giant DDRG, in which the newly restarted jets 
propagate through the bright cocoon formed by a previous 
active phase of its AGN. The current activity is indicated 
by the radio core. 

The maps of the outer and the inner radio lobes at five 
observing frequencies from 334 to 4860 MHz allow deter- 
mination of the spectral index distribution along the outer 
lobes, as well as extension of the spectral index frequency 
range in the inner lobes. The spectral index distribution in 
the outer lobes implies a distribution of the spectral (syn- 
chrotron) age of emitting particles along the ridge of these 
lobes, which increases from about 25 Myr at the edges of 
the extended old cocoon structure to about 65-75 Myr in 
the vicinity of the core. A velocity equivalent to the quo- 
tient of separation between the selected positions on the 
ridge and the difference in synchrotron age between these 
positions, which is higher than the advance speed of the 
jet's head, indicates a significant backflow in the old co- 
coon. Moreover, its average speed evidently decreases from 
the cocoon's edges towards the core. 

Using the DYNAGE algorithm (Machalski et al. 2007), 
we attempted to specify the dynamical age of both the outer 
and the inner pairs of radio lobes. Because each pair is 
characterized by a specific asymmetry in the lobes' length 
and brightness, we searched for a self-consistent solution of 
the analytical model, Demanding the same values for the 
jet power, central core density, and a comparable age for 
the opposite lobes, we find that 

(1) The age of the outer and the inner lobes is 132±28 
Myr and ~9±4 Myr, respectively. However the derived age 
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of the inner lobes is quite problematic because the almost 
pure power-law spectrum of its lobes prevents a unique age 
solution without any independent knowledge of either the 
restarted jet power or the ambient density within the old 
cocoon through which it propagates. This problem is solved 
with the assumption of limiting values for the latter param- 
eter discussed in detail in Sect. 5.3. 

(2) If the above assumption is correct, the jet powers 
during the initial and the restarted phase of activity are 
not similar. The restarted jet is one order less powerful 
than the original one. We think that this can support a hy- 
pothesis that the jet power is much more dependent (if not 
only) on the intermittent jet activity, than on the spin of 
a black hole. This intermittent activity can be connected 
with a stochastic transition between two accretion modes 
- the standard one where the angular momentum is trans- 
mitted outwards by viscous torques within a disk and - 
the 'magnetic' one, in which large-scale poloidal fields are 
developed (cf. Nipoti et al. 2005; Kording et al. 2006). 

(3) The magnetic field strength estimates along the 
ridge of the outer double structure vary slightly from the 
edges of the lobes towards the centre of the old cocoon. 
Unexpectedly, the field strengths in the inner lobes result- 
ing from the dynamical analysis seem to be twice lower 
than the corresponding values estimated in the spectral- 
ageing analysis. This discrepancy arises from the assumed 
different adiabatic indices in the equations of state: a non- 
relativistic one for the outer structure and the relativis- 
ts one for the inner double structure. The non-relativistic 
equation of state implies some energy dissipation process 
between magnetic field and thermal particles whose pres- 
ence in the old cocoon is expected. Therefore in the non- 
relativistic conditions for the magnetic energy evolution, a 
stronger magnetic field would be necessary to account for 
the observed brightness of the inner lobes. 
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Table A.l. Notations for the free parameters in the DYNAGE algorithm and other physical parameters used through 
the paper (excluding observational parameters of the radio source: D, Rt, V, P v , a u , defined in the text). 



Symbol 


Dimension 


Parameter 




Model 


free parameters to be assumed 


ao 


[kpc] 


radius of central core 


13 


[dim.less] 


exponent of ambient medium- density profile 




[dim.less] 


Lorentz factors of relativistic particles 


Tjct, Tb, r x , r c 

c 


[dim.less] 


adiabatic indices of the jet, magnetic field, ambient medium, and cocoon as a whole 


[dim.less] 


initial ratio of energy density of the magnetic field to that of particles 


k' 


[dim.less] 
[°] 


ratio of energy density of thermal particles to that of relativistic ones 


9 


orientation of the jet's axis 




Model 


parameters to be fitted 


Oinj 


[dim.less] 


injection spectral index 


t 


[Myr] 


dynamical age 


Qjet 


[W] 


jet power 


PO 


[kgm- 3 ] 


central core density 




Other 


physical parameters 


Psc.A, Psc.B 


[dim.less] 


exponents of ambient density profile in the self-consistent age solution A and B 


Phs, Pmin, Pc 


[Nm^ 2 ] 


hot-spot pressure, its minimum (equipartition) pressure, and cocoon pressure 


We, MB, Me 


[Jm- 3 ] 


energy density of relativistic particles, in magnetic field, and in cocoon as a whole 


Uout , U inn 


[J] 


total energy radiated from outer and inner lobes 


B, B e q v ,BiC 


[nT] 


strength of magnetic field, equipartition field, and inverse-Compton field 


k 


[dim.less] 
[m 3 ] 


ratio of energy density of relativistic particles to that of electrons 


rip, n c , n g 


proton, electron number density, and cold-gas number density 


m p 


[kg] 
[kgm 3 ] 


proton mass 


Pa 


ambient medium density 


(*I> 


[Myr] 


mean of t fit for the opposite lobes; i=>out, inn 


(Qjet.i) 


[W] 
[kgm- 3 ] 


mean of Qj ot fit for the opposite lobes 


(PO.i) 


mean of po fit for the opposite lobes 


(Pa,i> 


[kgm- 3 ] 


mean of p a derived for the opposite lobes 


r syn 


[Myr] 


spectral (synchrotron) age 




[GHz] 


frequency of spectral break 



